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Abstract 

In the case of the energy-momentum tensor related to "ordinary" matter (perfect fluid rep- 
resenting spin 1/2 and 1 fields), the equations of general relativity result in cosmological and 
gravitational collapse singularities — due to the fulfilment of both the strong and weak energy 
conditions. According to great unified theories, in superdense matter (both hot and cold), 
phase transitions take place, symmetry between the strong and electroweak interactions is re- 
stored/broken, and a scalar field is created/annihilated. In the scalar field regime, the strong 
energy condition is broken, but the weak one holds. However, the continuity condition for the 
pressure on the surface of a contracting star results in the occurrence of a compensating pseu- 
domatter field, for which both energy conditions are broken. On this basis, for a massive star 
predisposed to a gravitational collapse, a pulsation dynamics with no singularity is constructed. 
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Introduction 



In the case of the energy- momentum tensor related to "ordinary" matter (i.e., perfect 
fluid representing spin 1/2 and 1 fields), the equations of general relativity inevitably result in 
singularities — both cosmological (in the past) and gravitational collapse ones (in the future) 
[1-3]. The presence of singularities has provoked concern of physicists. Einstein himself had 
fought against them starting in 1918 [4]. 

The inevitability of singularities follows from the singularity theorems [3], which, in turn, are 
based on the strong and weak energy conditions. Denote by turn, A = 0, 1, 2, 3, the eigenvalues 
of the tensor T£, so that u>( ) = g (the energy density) and = —p(i), % = 1,2,3 (the 
principal pressures). The energy conditions read: g + YliP(i) — an d f? + > (strong); 
g > and g + P(i) > (weak). The strong energy condition relates to cosmological singularity, 
and both energy conditions are related to gravitational collapse. 

In the case of ordinary matter, g > and p^ > 0, so that both the strong and weak energy 
conditions are satisfied. 

The occurrence of singularities indicates a breakdown of general relativity: its dynamical 
equations are self-limiting in their scope [5] since they result in singularities, beyond which 
dynamics cannot be extended. 

According to grand unified theories, in superdense matter (both hot and cold), phase tran- 
sitions take place, after which symmetry between the strong and electroweak interactions is 
restored/broken and a scalar field ip ("extraordinary" matter) is created/annihilated [6]. In the 
scalar field regime, g = (l/2)y? 2 + V(tp), p w = p = (l/2)y? 2 - V(tp) [ = d/df, V{tp) is the 
potential], so that g > 0, g + P(%) > 0, g + YliP(i) = 2[y? 2 — ^(v 9 )]- Thus, the strong energy 
condition may be broken. Starting from this circumstance, in [7] a theory of deflation-inflation 
for the universe has been advanced, which results in the contracting-expanding universe with 
no singularity. 

On the other hand, for the scalar field, the weak energy condition holds, so that it seems that 
gravitational collapse (i.e., singularity) is inevitable. However, the existence of the boundary 
surface of a contracting star results in the occurrence of an external pseudomatter field [8], 
which compensates for the discontinuity of the scalar field pressure on the boundary surface. 
Thus, extraordinary matter is comprised of the interior scalar and exterior pseudomatter fields. 
(The latter represents dark matter.) 

For the pseudomatter field, the result is this: P( n0 rmai) < 0, g < 0, so that both energy 
conditions are broken, which makes it possible to avoid singularity. 

In the present paper, for a massive star predisposed to a gravitational collapse, a pulsation 
dynamics with no singularity is constructed. 
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1 Basic equations 

1.1 The extended Einstein equations. Matter and pseudomatter 

In the extended Einstein equations [8] 



1.2 Spherically symmetric star. Dimensionless quantities 

We consider a spherically symmetric star; denote by a its radius and introduce dimensionless 
quantities (the subscript dim stands for dimensional): 

time: t = £dim/ a ; radial coordinate: r = r d [ m /a; 

radius of sphere with radial coordinate r: K(r,t) = K dim (r,t)/a; 

metric: ds 2 = ds^ im /a 2 ; 

the Einstein tensor: G 1 ^ = a 2 G^ dim ; the energy-momentum tensor: = a 2 xT^ dim ; 
energy density: g = a 2 xpdi m ; pressure: p = a 2 xpdi m ; 
scalar field: tp = V a 2 >np dim ; potential: V = a 2 xVdi m - 

1.3 Metric, the Einstein equations and tensor components 

In the case of spherical symmetry, in the synchronous reference frame, metric is of the form [2] : 




V 



(1.1.2) 



T ps ^ compensates for the discontinuity of T m ^ in time and space. 



ds 2 = dt 2 - e A ^dr 2 - K 2 (r, t)(d6 2 + sin 2 9 d</> 2 ) 



(1.3.1) 



Again, in the case of spherical symmetry, the nonvanishing tensor components are 



A 0> A l> A 2 — A 3> A 0> A l> A — ^5 J m, Jps 



(1.3.2) 



where 



x° = t, x 1 = r, a; 2 = 6 
Thus, (1.1.1) in the dimensionless form reduces to 



9, x 3 = (j) 



(1.3.3) 



G\ - 8irTl = 



(1.3.4) 



G\ - 8nT 2 = 
G° Q - 8ttT ° = 
G\ - 8irT? = 



(1.3.5) 
(1.3.6) 
(1.3.7) 
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Now [2], (• = d/dt, 1 = d/dr) 



~i 2K K 2 1 (K') 2 , 
3*4 + * + ™ (1 .3.9) 



AT 2 4 2K 2 \ K K 

1.4 Matter and pseudomatter equations 

Matter and pseudomatter equations are of the form 

T^ ;A = 0, z = m,ps (1.4.1) 

i.e., 

K x + rlT z ;-r; x T z x a = o (1.4.2) 

The nonvanishing Christoffel symbols are [2] 



r 1 
1 11 


= A'/A 


r 1 - 

1 22 _ 


: —e~ x KK' Tl 3 = —e~ x KK' sin 2 9 


r} = a/2 


r 2 

1 12 


= K'/K 


r 2 

1 20 


= K/K r 2 3 = -sin^cos^ 




r 3 

1 13 


= K'/K 


r 3 

1 30 


= X/X r 3 3 = cot e 




r° 
1 ii 


= Ae A /2 


r° 

1 22 


= KK r° 3 = KK sin 2 9 





(1.4.3) 

As a consequence of (1.4.1), the dynamical equations are (1.3.8), (1.3.9), whereas (1.3.10), 
'1.3.11) are constraints on initial conditions. 



2 Setting of a problem 

2.1 Metric singularity 

The definition of the singularity of metric is this: 

K(r, t) -> and e A(r '* } ^0 for t -> t sing < oo (2.1.1) 

The existence of singularities is the consequence of the fulfilment of the energy conditions [3] . 

A gravitational collapse results in the singularity (2.1.1). There arises the problem of a mass 
point with a metric singularity at the point (rather than a mass point in a regular metric field). 
For this problem, the equations of general relativity, specifically the geodesic equation, fail. 
Thus, it is impossible to construct a consistent dynamics for the universe with singularities. 

In this connection, we quote Penrose [5]: 
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"This is not to say that some mathematically precise concept of 'singularity' should now 
form part of our description of physical geometry — though much elegant work has been done 
in this direction in recent years. Rather, it seems to be that the very notion of spacetime 
geometry, and consequently the physical laws as we presently understand them, are limited in 
their scope. Indeed, these laws are even self- limiting, as the singularity theorems ... seem to 
show. ...There is a need for new laws in any case..." 

2.2 The energy conditions 

Denote by w the eigenvalues of the energy-momentum tensor, i.e., solutions to the equations 

T;v v = wv^ w = w (x) , A = 0,1, 2, 3, v^O (2.2.1) 
The standard notation is this: 

W(P) = Q, W(i) = ~P(i), « = 1,2,3 (2.2.2) 

The eigenvalue g may be interpreted as the density, and the eigenvalues are called the 
principal pressures [3]. 

The strong energy condition is 

3 

Q + J2PW-° and £ + P(i)>0, i = 1,2,3 (2.2.3) 
i=i 

The weak energy condition is 

g>0 and g + p {i )>0, i = 1,2,3 (2.2.4) 

The fulfilment of the strong energy condition results in a cosmological singularity. The 
fulfilment of the weak or strong energy condition results in a gravitational collapse singularity 
[3]. 

The case of cosmology has been considered in [7]. In the present paper, we are interested 
in the case of a star. Therefore our problem is, first of all, to determine if both the weak and 
strong energy conditions may be broken. 

2.3 Ordinary and extraordinary matter 

With respect to density and pressure, there are ordinary matter (perfect fluid representing spin 
1/2 and 1 fields) and two kinds of extraordinary matter: scalar field and pseudo matter. 

Ordinary matter 

In the case of ordinary matter, 

£?ord>0, p ord >0, 2 = 1,2,3 (2.3.1) 
so that both the strong and weak energy conditions are fulfilled. 
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Scalar field 

According to grand unified theories, in superdense matter (both hot and cold), phase tran- 
sitions occur, after which symmetry between the strong and electroweak interactions is re- 
stored/broken, with the result that a scalar field ip (extraordinary matter) is created/annihilated 
[6]. The energy density and pressure of a spatially homogeneous scalar field ip(t) with potential 
V((p) are of the form 

1 1 

Qsc = -0 2 + V((p), Psc(i) =p sc = -0 2 - V(<p), i = 1,2,3 (2.3.2) 

We have 

Qsc > (2.3.3) 
Qsc + Psc^) = V? 2 > 0, i = 1,2,3 (2.3.4) 
Qsc + ^Pscii) = Qsc + 3p sc = 2[y? 2 - V(tp)\ (2.3.5) 

i 

Thus, the strong energy condition may be broken, but the weak energy condition is fulfilled. 

Pseudomatter 

It is the behavior of the energy conditions for pseudomatter that is the crux of the singularity 
problem. 

2.4 The problem 

Now the problem is this: To determine if due to the involvement of pseudomatter both the 
strong and weak energy conditions may be broken and, if so, to construct a pulsation dynamics 
of a star. 

3 Model 

3.1 Interior and exterior 

Assume that in the interior of the star there is only matter and in the exterior only pseudo- 
matter: 

interior: < r < 1, T ps » = 0, T; = (3.1.1) 
exterior: 1 < r < oo, T m ^ = 0, T; = T ps » (3.1.2) 

3.2 The matter energy- momentum tensor and matter equations 

Let 

^mo =: Qnn ^ml = ^~m2 = ^~m3 =: Pm (3.2.1) 

and 

T m ? = Tj = (3.2.2) 

This corresponds to the perfect fluid: 

T m ^ = (ftn + Pm)Mm /1 M m i/ - pS^ (3.2.3) 



with 

Um 1 = 0, tt m0 M m ° = 1 

Now matter equations (Subsection 1.4) reduce to 

T m 2^x = T m3 f x = 

T mQ + (r 10 + 2r 20 )(T m0 — T ml ) = o 

T m i]i — 

Thus the pressure is spatially homogeneous: 

Pm = 0, Pm = Pm(t) 

3.3 Homogeneous density, interior metric, 

the interior Einstein equations, and matter equation 

In view of (3.2.8), put 

£m = 0; ftn = ftn(t) 

Under the homogeneity conditions (3.3.1), (3.2.8), we have in the interior 



K int (r,t) = R(t)r, (e A ) int = 



1 _ r 2 



, < r < 1 



and 



(ds 2 ) int = dt 2 - R 2 {t) 



dr 2 



I _ r 2 

Now the interior Einstein equations reduce to 



+ r 2 (d# 2 + sin 2 fld0 2 ) 



< r < 1 



2R R 2 1 

3 1 W + TP 1 " a7II!m = 



and the matter equations reduce to 



QmR + 3(g m +p m )R = 



(3.2.4) 

(3.2.5) 

(3.2.6) 
(3.2.7) 

(3.2.8) 



(3.3.1) 

(3.3.2) 
(3.3.3) 

(3.3.4) 
(3.3.5) 



(3.3.6) 



Dynamical equations are (3.3.4), (3.3.6), and (3.3.5) is a constraint on initial conditions. 
Without pseudomatter we would have 



which would result in 
Later on, we will put 



Pm(M) =0 

Pm(r,t) = 

Pord = 



(3.3.7) 

(3.3.8) 
(3.3.9) 
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3.4 Relation between matter density and pressure 

For ordinary matter, the relation between density and pressure is given by the state equation: 

Pord = Pord(ford) (3.4.1) 

For scalar field, density and pressure are related via the field: 

Q B c = \<i> 2 + V(<p), p sc = ^ 2 -%) (3.4.2) 
The matter equation (3.3.6) reduces to 

rlT/1 

= (3.4.3) 



i.e., either 



dV 

R(p + 3R(p + R— 
dip 



ip = (3.4.4) 



or 

d^ 

R(p + 3R(p + R— = (3.4.5) 

dip 

If both (3.4.5) and (3.4.4) hold, then 

d^ 

- = (3.4.6) 

In the case of (p — 0, 

- P sc = Q sc = V(<f) = const (3.4.7) 

3.5 The pseudomatter energy- momentum tensor, 

pseudomatter equations, and the exterior Einstein equations 

Taking into account matching conditions at r = 1 [2] and (3.2.1), (3.2.2), (3.2.3), (1.4.1), we 
put 

Tps2 = ^ps3 = 0, T psl = Tp SQ = (3.5.1) 

but must retain 

Tpsi = - Ppsi? -^pso =: Qps (3.5.2) 
The pseudomatter equations take the form 

^pso°o + ^ + 2— ^ Tp S o - = (3.5.3) 

T psl \ + 2^T ps \ = (3.5.4) 
and the exterior Einstein equations are these: 
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G\ - 8nT ps \ = (3.5.5) 

G 2 2 = (3.5.6) 

G° - 8nT ps ° = (3.5.7) 

G° = (3.5.8) 

with G£ given by (1.3.8)-(1.3.11). 

(3.5.3)-(3.5.6) are dynamical equations; (3.5.7), (3.5.8) are constraints on initial conditions. 

3.6 Matching conditions at r = 1 

Matching conditions at r = 1 are these [2]: 

Ppsl (l,t)=p m (l,t) (3.6.1) 

K ext (l,t) = K int (l,t) = R(t) (3.6.2) 
From (3.6.1), (3.6.2), (1.3.4), (1.3.8) follows 

[(K') 2 e- X } cxt (l,t) = [(Jf) s e-*WM) (3.6.3) 

Now, 

[(K') 2 e- X ]' = K'[(2K" - K'\')e- X ] (3.6.4) 

so that the expression (2K" — K'X) involved in (1.3.9), (1.3.10) does not contain 5(r — 1). 
Next, 

[(K') 2 e- X \ = K'e~ x [2k' - \K'\ (3.6.5) 
which, with (1.3.11), (3.5.8), results in 

[{K') 2 e- X \ = (3.6.6) 

Again, 

l(K>) 2 e-%t = 1 - r 2 , [(K') 2 e-% t (h = (3.6.7) 

Thus, 

[(K>) 2 e- x Ua,t)=0 (3.6.8) 

so that 

iCt(M) = (3.6.9) 
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3.7 Pseudomatter pressure and density 

The equation for p psl is (3.5.4): 



K' 

Ppsi + 2 

whence 



+ 2— P P si = (3.7.1) 



Ppsl K 2 = /(t), Ppsl (r,t) = -^j^ (3.7.2) 
From (3.6.1) follows 

jMd =p - m (3 - 7 - 3) 

Thus, in view of (3.6.2), 

ft.i(r,f) = p m (M) J^Tjy = M*)j^j (3-7-4) 

Note: as long as 

lim tf(r,t) = (3.7.5) 

(Ppsi)int is inadmissible. So 

P P si(r,0-Pm(O x f ^,t) (3 ' 7 - 6) 
Next, the equation for g ps is (3.5.3): 

Q ps + ^ + 2^j £? ps + ^p ps i = (3.7.7) 

Denote by ib-r the instant of the symmetry broken-restored transition and by t r _b the instant 
of the reverse transition (if the latter occurs, i.e., if there is no singularity). Put 

^ps(ib-r) = (3.7.8) 

then 

1 f de x(r ' t)/2 

MM) = - x2(M)eA(r , t)/2 J -^^ Pm (t)R 2 (t)dt (3.7.9) 

*b-r 

3.8 Breakdown of the energy conditions 

Assume that there is no singularity. Then for 

*b-r < t < t z - b (3.8.1) 

we have 

Pm = Psc = -<P 2 - V(<p) (3.8.2) 
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Let 

then, by (3.7.6), 
Introduce 

We have 
Thus 



p sc < (3.8.3) 
Ppsi < (3.8.4) 



^b— r ~t~ ^rb 



(3.8.5) 



Q e \/2 Q e X/2 

— (r,* )=0, — ( r ,t)^0 for t$ t (3.8.6) 



£ ps < for t < t r _ b (3.8.7) 
We have obtained the following results: for t given by (3.8.1) 

£ps + ^Ppsi < and ^ P s + Ppsi < (3.8.8) 

i 

g ps < and g ps + p psl < (3.8.9) 

so that both the strong and weak energy conditions are broken. 
In the interior, the strong energy condition is broken as long as 

£? sc + j^p^ = &c + 3 ^ c = 2 & 2 - v ^ < ( 3 - 8 - 10 ) 

i 

For the metric (3.3.3), it is the breakdown only of the strong energy condition that is essential. 

Thus, the energy conditions are broken in the whole spacetime, and our assumption is 
justified. 

Now we may construct a pulsation dynamics, i.e., one free of singularity. 

4 Pulsation dynamics 

4.1 Broken symmetry (ordinary matter) stage 

Consider a cycle of a pulsation dynamics. Put 

— initial = tfinal > (4-1.1) 

in which case 

- t b _ r = t r _ b > 0, t = * b " r | tr ' b = (4.1.2) 
The broken symmetry, i.e., ordinary matter stage takes place for 

tE [initial, tb-r) ^J(tr-b, tfinal] (4.1.3) 

with 

^initial) = ^(ifinal) =0, X = K,\ (4.1.4) 
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Interior 

The dynamical equations are: 



with 

and the constraint 



2R R 2 1 

Q OTd R + 3(£Wd + Pord)R = 
Pord = Pord(^ord) 



R 2 1 
3 ' R 2 + R 2 ' ~ 87rf?ord = ° 



Exterior 

The dynamical equations are: 

^ + 8vrp ord (t)j|^ = 
G* = 

and the constraints are: 

Go " 87r P P s = 
G° 1 = 

[G^ are given by (1.3.8)-(1.3.11)]. 
In accordance with (3.7.8), we put 

Pord = 0, p ps = 0, g ps = 

In this case, there exists an exact, well known solution to the above equations [2]. 
Matching 

Matching conditions at r = 1 are determined in Subsection 3.6. 

4.2 Restored symmetry (scalar field) stage 

The restored symmetry, i.e., scalar field and pseudomatter stage takes place for 

t E [tb- r , t r _b] 

Interior 

The dynamical equations are: 



2R R 2 1 



\v 2 ~ V( V ) 



dV 

<p = or R(p + 3R(p + R— = 

dtp 
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with the constraint 



, R 2 1 , 

3| ^ + ^' " " r 



if + V(<p) 



From (4.2.2), (4.2.4) follows 
where 

The inequality 

must hold. 
We have 

and 



r - t 2 r = o 



r 2 = f [V( V )-^] 



T 2 > for * (4.2.1) 



R(0) = 0, 0(0) = 



R(0) = 



r(0) V 8nV(ip(0)) 



Exterior 

The dynamical equations are: 



-f - V{tp) 



Gj = 



and the constraints are: 



Gq - 8ng ps 
G? = 

with G^ given by (1.3.8)-(1.3.11) and g ps by (3.7.9). 

Matching 

Matching conditions at r = 1 are determined in Subsection 3.6. 

4.3 Matching conditions at phase transitions 

Matching conditions at the phase transitions (t = tb-r,^r-b) are these: 

f?ord(£b-r/r-b) = f?sc (^b-r/r-b) 

i^b-r/r-b + 0) = i2(t b -r/r-b - 0) 
R(h-r/r-b + 0) = i?(t b _ r / r _ b — 0) 



(4.2.4) 

(4.2.5) 
(4.2.6) 
(4.2.7) 

(4.2.8) 
(4.2.9) 



(4.2.10) 

(4.2.11) 

(4.2.12) 
(4.2.13) 



(4.3.1) 

(4.3.2) 
(4.3.3) 
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4.4 The simplest model: Constant scalar field 

Consider the simplest model: 



in the restored symmetry stage. 
We have 



with 



and the initial conditions 



The solution is this: 



or, in the dimensional quantities, 



(p = 0, <f — const 



R - TR = 



87rV((p) 



const 



R(0) 



r' 



R(t) 



R 



dim 



R(0) = 



coshTt 



cosh r dim t di: 



dim 



/ 87rxV(y) 

1 dim — V o , ^dim 



V 



4.5 Pulsation dynamics stages (one cycle) 

The pulsation dynamics stages for one cycle are these: 

initial : R = -Rmax, R = 0, R < 
(initial, ^b-r) : R < 0, R < 

(t b _ r ,0): i?<0,i?>0 

t = 0: R = R min ,R = 0,R> 

(0,* r _ b ): i?>0,i?>0 

(*r-b, ^final) = -R > 0, i? < 
£final : -R = -Rmax, R = 0, R < 



(4.4.1) 

(4.4.2) 
(4.4.3) 
(4.4.4) 
(4.4.5) 

(4.4.6) 
(4.4.7) 



Acknowledgments 



I would like to thank Alex A. Lisyansky for support and Stefan V. Mashkevich for helpful 
discussions. 



14 



References 

[1] L.D. Landau, E.M. Lifshitz, The Classical Theory of Fields (Pergamon Press, 1975). 

[2] Hans Stephani, Relativity (Cambridge University Press, 2004). 

[3] Robert M. Wald, General Relativity (The University of Chicago Press, 1984). 

[4] A. Einstein, Kritisches zu einer von Hrn. de Sitter gegebenen Losung der Gravitationsgle- 
ichungen (Sitz. Preuss. Akad. Wiss., 1918, 1, 270-272, 1918). 

[5] R. Penrose, Singularities and time-asymmetry, in: General relativity, An Einstein cente- 
nary survey, ed.: S.W. Hawking, W. Israel (Cambridge University Press, 1979). 

[6] A.D. Linde, Particle Physics and Inflationary Cosmology (CRC Press, 1990), ("Nauka", 
Moscow, in Russian, 1990). 

[7] Vladimir S. Mashkevich, The Eternal Closed Universe: Deflation-Inflation (arXiv: 
0910.0894v2 [physics. gen-ph] , 2009). 

[8] Vladimir S. Mashkevich, On Dark Matter Problem: Pseudomatter — Concept and Appli- 
cations, arXiv: gr-qc/0902.3508 (2009). 



15 



